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Abstract 

We consider the M — > limit for tagged particle diffusion in a 1-dimensional 
Rayleigh-gas, studied originaly by Sinai and Soloveichik ^3], respectively by Szasz 
and Toth [Ej. In this limit we derive a new type of model for tagged paricle diffu- 
sion, with Calogero-Moser-Sutherland (i.e. inverse quadratic) interaction potential 
between the two central particles. Computer simulations on this new model repro- 
duce exactly the numerical value of the limiting variance obtained by Boldrighini, 
Frigio and Tognetti in [3J. 

1 Introduction 

The problem of deriving relevant information on the diffusive scaling limit of tagged 
particle motion (i.e. self-diffusion), from microscopic principles has been undoubtedly at 
the heart of mathematically rigorous statistical physics of time dependent phenomena, 
at least since Einstein's groundbreaking work. Mathematically rigorous investigation of 
tagged particle diffusion in systems of particles governed by deterministic (Hamiltonian) 
dynamics is notoriously difficult even in one dimensional models. After remarkable ad- 
vances achieved up to the late 1980-s (see Section 01 below and references cited there), in 
the last twenty years there seemed to be less intense activity in the field. This is certainly 
due to the difficulty of these problems and lack of technical tools to attack them. 
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In the present note we make a small but hopefully not completely irrelevant contri- 
bution to the subject. We investigate the M — > small mass limit of the tagged particle 
diffusion in the so-called 1-dimensional Rayleigh-gas. This system consists of an infinitely 
extended one-dimensional system of point-like particles of mass 1 and one single tagged 
particle of mass M immersed in it. The particles perform uniform motion and interact 
through elastic collisions. The system is distributed according to the equilibrium Gibbs 
measure. This means independent exponentially distributed inter-particle distances and 
independent normally distributed velocities with mean zero and inverse mass variances 
(that is: equidistributed kinetic energy). It is a fact that the dynamics of the infinitely 
extended system is almost surely well defined under this stationary measure. That is: no 
multiple collisions and no accumulation of infinitely many particles in finite time occurs. 
Randomness comes into the problem only through the thermal equilibrium of the initial 
condition, otherwise the dynamical evolution is deterministic. The central question is un- 
derstanding the diffusive scaling limit of the trajectory of the tagged particle: A~ x l 2 QAt-, 
as A —>■ oo. There exist a number of deep and interesting results related to this problem 
which will be shortly surveyed in Sectional In the present note we investigate the limit 
when M < 1. We prove that in this limit the system becomes equivalent to another, 
new model of tagged particle motion, differing from the one described above in having 
instead of one central particle of different mass, all particles of the same mass but the 
two central particles interacting via a Calogero-Moser-Sutherland-type repulsive poten- 
tial, with random strength parameter. This result explains some phenomena observed in 
earlier computer simulations on the Rayleigh-gas. In particular the instability observed 
for small values of M. We also present numerical simulations on this new model. Our 
simulation results reproduce very accurately the numerical value of the limiting variance 
of tagged particle in the Rayleigh-gas, in the M — > limit, which were obtained in 
We claim that our result not only reproduces (from a completely different approach) the 
numerical value but also gives theoretical explanation of the phenomenon. 

The paper is organised as follows: In Section |21 we define the models of interacting 
particle systems considered, their stationary Gibbs measures and the stochastic processes 
whose diffusive asymptotics is later analysed. In Section El we briefly survey the existing 
earlier results (rigorously proved and numerical, alike) on tagged particle diffusion in 1-d 
Rayleigh-gas. In Section HI we properly state and prove the theorem which states that 
in the M — > limit the 1-d dynamics of the Rayleigh-gas with tagged particle of mass 
M converges (trajectory-wise, in a natural topology) to the dynamics of the 1-d gas of 
particles with Calogero-Moser-Sutherland interaction between the two central particles. 
Finally, in Section we present our new numerical results referring to this new type of 
interacting particle system. We should emphasize here that our numerical results are 
not just accurate reproduction of older computer experiments but are performed on a 
genuinely new type of model. One of the main points of this paper is exactly the fact 
that these genuinely new numerical results are in accurate agreement with the results of 
Boldrighini, Frigio, Tognetti giving independent enhancement and theoretical expla- 
nation to them. 
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2 Models: state space, dynamics, stationary mea- 
sures 



In this section we describe the models considered throughout the paper. In section I2~T1 
we present a formal definition of the state spaces and the natural measures on them. 
Section EIH which gives a more verbal description of the time evolution in our dynami- 
cal systems, clarifies that these models indeed correspond to the one dimensional gases 
mentioned in the Introduction. 



2.1 State spaces and stationary Gibbs measures 

Let 

Q ± := {cj ± = (x±i, v ±i )Zi ■ (x±i, v±i) 6^x1, x±i = 0, ±(x ±(i+1) - x ±i ) > 0}. 

With slight abuse of notation and terminology sometimes we don't distinguish between 
uj ± and the (unordered) set of points {(x±i, v±i) : % — 1, 2, . . . } C R ± x R. We endow the 
spaces fi ± with the topology defined by pointwise convergence: lu^ — > u; ± iff (x±i, v±i) n — > 
(x±i, v±i), for all i — 1, 2, . . . . This is a metrizable topology and makes complete and 
separable (i.e. Polish) spaces. 

We denote by fi^ the following probability measures over fl^, respectively. Under ^ 
the random variables := ±(x±(j + i) — x±i), r]±j := v±j, i,j = 1, 2, ... , are completely 
independent, with exponential, respectively, normal distributions: 



P(C±i e (x,x+ dx)) = l {x > 0} e x dx, P(r]±j e (v,v+ dv)) = 




We shall consider two different types of particle systems in this paper. Their state 
spaces will be 

Q 1 ■= {(uj + ,uj-,z,u,V) : G Q ± , zeR+, ue[-l,l], V G R}, 
Q 11 ■= {(uj + ,uj-,z) : u± G Q ± , zeR+}. 

We also define the natural projection between these spaces: 

II : Q 1 -> Q 11 , n(cu + , uj~, z, u, V) := u~,z). 

In order to define the relevant probability measures on the state spaces Q 1 and Q 11 
first we introduce some notation. Let the random variables W and ( be independent and 
distributed as a standard Gaussian, respectively, as a standard T(2). Let 72(2) be the 
density of the distribution of (, g(c) the density of the distribution of \W(\, and <p c ( z ) 
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the density of the conditional distribution of (, given \W(\ = c: 



72(2) := ze 
g(c) 



^ c 2 
-exp{-^ - — 
g(c) V 7r 2z A 



l2{z) = / (p c {z)g(c)dc. 



Clearly, 



The probability measures considered on the state spaces Q 1 and Q 11 are fi I,M (defined 
on Q 1 ) which depends on the positive parameter M, respctively, fi II,c (defined on fl n ) 
which depends on the positive parameter c, and finally fi 11 (also defined on Q 11 ) which is 
a mixture of the measures [i n,c : 



//' "(do/) := x //-(do;") x l2 (z) dz x - du x W ^ e w2/2 dV, (1) 

2 V 27r 



t J ' M (do/) 



/x Ji ' c (du/ J ) := /i + (dcj + ) x jr(cLr) x v? c (z)dz, (2) 
H n {duj n ) := /i+(do; + ) x 0-(daT) x 7a (z) dz 

fjc n > e {du 1I )Q(c)dc. (3) 



The measures // ,M , respectively, // /,c will be the natural Gibbs measures corresponding 
to the dynamics of our systems, to be defined in the next subsection. 



2.2 Dynamics 

We define the dynamics of the systems considered verbally, rather than writing formulas. 
The two types of dynamics considered will be called of type I, respectively, of type II. 
Their state spaces will be Q 1 , respectively, Q n . These will actually be families of dynamics 
parametrized by the fixed parameters M > 0, respectively, c > 0. 



2.2.1 Dynamics of type /: 

For precise formal definitions and basic facts about these dynamics see jHj, 120] ■ The 
system consists of particles indexed . . . , —2, —1, 0, +1, +2, .... The system is observed 
from the tagged particle of index 0. The tagged particle has mass M, the other particles 
have unit mass. Positions and velocities of the particles in the system are encoded in 
(u + , u~ , z, u, V) as follows: V is the velocity of the tagged particle, x±i ± z(l ± u)/2 and 
v±i is the position, respectively, the velocity of the particle of index ±i, i — 1,2, ... . The 
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untagged gas particles perform uniform motion on the line and don't interact between 
themselves, when two of them meet and cross each other's trajectory they exchange their 
index. The tagged particle, while isolated from the others, also performs uniform motion 
and collides elastically at encounters with an untagged gas particle. At these collisions 
the outgoing velocities v ont ,v out are determined by the incoming velocities V m , v 111 as 
follows: 

M - 1 • 2 f 2M ■ M - 1 • 

V out = V m + v m v ont = V m v m (A) 

M+l M+l ' M+l M+l K J 

Mind that the untagged gas particles never exchange their order with the tagged particle 
and the index +i of a particle denotes its actual relative order with respect to the tagged 
particle. 

The measure fj, ' defined in is Gibbs measure for this dynamics, invariant for 
the system as seen from the tagged particle. It is a fact (see [Hj, [TH], (201) that the 
dynamics is /i 7 ' M -a.s. well defined: starting the system distributed according to /i /,M , 
with probability 1 no multiple collisions will occur and the system remains locally finite 
indefinitely. We denote by Sl' M the measure preserving flow defined by this dynamics 
on (£y,//> M ). 

The velocity and displacement process of the tagged particle is 

V/'» = V^iui 1 ) := VWcS) 
Q\> M = Q[ M ^):= /V^(u/)ds. 



In section 0] it will be more convenient to describe the dynamics from a fixed exterior 
point of observation. The absolute locations of the gas particles in the system as seen 
from such a fixed exterior frame of reference are 

V±i{t) ■= Q[' M + x ±i (t), where x ±i {t) := x ±i {Sl ,M u 1 ), i = 1, 2, . . . . (5) 
We also introduce the variables 

V^ M = ^(u/) := \ (y^S^+v^S^u; 1 ) 

These are the velocity and position processes of the centre of mass of the particles next 
to the right and to the left of the tagged particle. We need the position process Q\' M for 
later comparison with a similar process defined for the dynamics of type il in the next 
paragraph. Mind that the random process t i— > (Q t ,M — Q\ ,m ) is stationary and thus tight, 
uniformly for t > 0. As a consequence, (Ql' M / Vi — Q{ ,M / \^t) — > in /i 7,M -probability 
(actually //' M -a.s.) as t — >• oo. 
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2.2.2 Dynamics of type 77: 



The system consists of particles of unit mass indexed . . . , — 2, — 1, +1, +2, . . . . Mind 
that there is no particle of index in this system. The system is observed from the 
centre of mass of particles of index +1 and — 1, we call this the central observation 
point. Positions and velocities of the particles in the system are encoded in (u> + , uj~, z) 
as follows: x±i ± z/2, respectively, v±i is the position relative to the central observation 
point, respectively, the velocity of the particle of index ±i, % = 1,2,.... Clearly, z 
denotes the distance between the two central particles of index +1 and —1. Particles 
move uniformly on the line except for the two central particles of index +1 and —1 which 
interact via the inverse quadratic pair potential U(z), or equivalently repelling force F(z): 



where c 2 > is a fixed parameter and z is the distance between the two central particles. 
When two gas particles meet and cross each other's trajectory they exchange their index. 
But mind that due to the strongly repulsive interaction between the two central particles, 
these two will never meet and thus particles will never change the sign of their index. 
The index ±z of a particle denotes its actual relative order with respect to the central 
observation point. 

Remark. In the literature of completely integrable Hamiltonian systems the pair po- 
tential © is usually called Calogero-Moser-Sutherland interaction and leads to one of 
the most notorious completely integrable 1-d systems, see [Sj, [H] and [THJ for the first 
original publications. 

The measure fi II,c defined in (j2J) is Gibbs measure for this dynamics, invariant for the 
system as seen from the centre of mass of the two central particles. It is again a fact that 
this dynamics is // /,c -a.s. well defined. We denote by iS/ 7 ' c the measure preserving flow 
defined by this dynamics on (Q 11 , fi 11 ,c ) . 

The velocity and displacement process of the point of observation is 



Again, absolute locations of the gas particles in the system as seen from a fixed exterior 
frame of reference are expressed similarly to (JHJ). 




(6) 



V t ' = V t • (u ) := ' w )+v +1 (S t uj 



)) 
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2.2.3 Stochastic processes considered 

In this paper we consider the following stochastic processes: 



Qt' M 


= Qt ,M 




with random u 1 distributed according to 




Qt' M 


= Q\ M 


OA 


with random u 1 distributed according to 




Q"' c 


= Q"' c { 


<A 


with random to 11 distributed according to 




Q" 


= Qf' c ( 




with random (uj 11 , c) distributed according 


; to dfi II,c g(c) dc 




= Qf' c , 


with random c distributed according to g(c) dc. 



This means that the process Qf is a g(c) dc-mixture of the processes Q t 



3 Survey of earlier results 

In this section we summarize the old results - rigorously proved and numerical - regarding 
various limits for the motion of the tagged particle in the model of type /. In Section H] 
we formulate and prove a new result concerning the M — > behaviour of these systems. 
In Section |3] we describe our new numerical results, referring to this new model. 

In all cases we are interested in the diffusive scaling limit of the displacement of the 
tagged particle motion, that is in the asymptotics of the rescaled process 

t ^ A~ 1/2 Q At , as A -> oo. 

Throughout the paper we denote by A this scaling parameter. 

We shortly survey the existing results on the asymptotics of the tagged particle motion 
in model of type / in historical order. The constants 

a 2 := 0.627..., a 2 := w 0.798 . . . 

will play a key role in the formulation of these results. 

3.1 The M = 1 case: 

The case when the tagged particle has the same mass as the rest of the gas particles was 
investigated and solved in Spitzer (1969) [T7j. For the roots of these ideas see also Harris 
(1965) jH]. In P?] the following invariance principle is proved: 

for M = 1 : A- 1/2 Q^ aW t , as A -> oo, 

where ^> stands for weak convergence of the sequence of processes (see [2J for weak 
convergence of processes), and Wt is a standard 1-d Borwnian motion. That is: WWt is 
a Brownian motion of variance a 2 . 
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3.2 The Ornstein-Uhlenbeck limit: 

Holley (1971) [0] considers the following limit when the mass of the tagged particle is 
rescaled in the same order as the time scale factor. Let m e (0, oo) be fixed. Then 

for M = mA: {A 1 ' Vjf , A'^Q 1 ^) =► (t£\ C) , as A - oo, 

where J)™ and £|™ are the Ornstein-Uhlenbeck velocity, respectively, position processes 
defined by the SDEs 

d V ? = - 7 (m)C dt + dW t , C = / V?ds, 

Jo 

with friction and dispertion parameters 

j(m) ■= — J-, D{m) := -^J -. 

For a version in higher dimensions of this type of result see Diirr, Goldstein, Lebowitz 
(1981) 

It is important to remark (sec |20j), that 

£™ aW t , as m -> 0. 

This means that taking first Holley's limit, then m —>■ we obtain a Wiener process of 
variance a 2 as the diffusive scaling limit of the displacement of the tagged particle. 

3.3 Bounds for the limiting variance for any M: 

Sinai, Soloveichik (1986) jEj, respectively, Szasz, Toth (1986) [TH] consider the case of 
arbitrary fixed mass M of the tagged particle. In these papers very similar rersults are 
proved in completely different ways. The results are summarized as follows: 

for M < A : ah < lim inf Var (A _1/2 Q^ M ) < lim sup Var (A~ 1/2 Q^ M ) < a 2 t. 

Mind that these bounds are independent of the mass of the tagged particle. For surveys 
of these results see also [T3], [T3], [T£l, pi] . 

Any rigorous result regarding the mass dependence of the limiting variance 

ai f :=hmVar(r^Q^) 

t— >oo 

remains one of the most interesting open questions in this context till today. The only 
known case is Spitzer's result <j\ = a. For numerical results see subsection IH.5I below. 
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3.4 Large mass Wiener limit: 



In order to interpolate between the M = const, cases (see subsection 13.3)1 and Holley's 
limit (see subsection 13.2) ) Szasz and Toth (1987) j2D] considered the limit with asymp- 
totics 1 <C M <C A, as A —>■ oo. Here the main result is the following invariance principle: 

for A 1/2+e < M < A : A'^Q 1 ^ aW t , as A -> oo. (7) 

Actually, the scaling limit (j2J) should hold for 1 C M C A but the method of proof in 
|2U] based on a coupling argument breaks down for 1«M< A l l 2+£ . For a survey of 
the results recalled in this and the previous paragraph see also fU\ . 



3.5 Earlier numerical results 

Following jHj and ^J] various numerical investigations were performed in order to es- 
tablish the mass dependence of the limiting variance: M > a 2 M . 

The relevant numerical investigations performed in the late eighties, early nineties 
are published in Omerti, Ronchetti, Diirr (1986) ^2], Khazin (1987) ^0], Boldrighini, 
Cosimi, Frigio (1990) jB], Fernandez, Marro (1993) 7J. These results clearly suggest the 
qualitative dependence M \— > a 2 M shown in Figure [T] 




Figure 1: Qualitative dependence M \— > a 2 M suggested by earlier numerical works 

As about the M — > limit: in all these papers it is remarked that the numerical 
simulations for small mass of the tagged particle are unreliable due to instability. On 
the other hand there was agreement between all researchers intrested in these questions 
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that \im M ^ a 2 M = a 2 ~ 0.798 . . . should hold. The "starightforward argument" was the 
following: the tagged particle of extremely small mass must have very small effect on the 
system, vanishing as M — > 0. So, in the M — > limit the displacement of any marked 
particle (in particular the one next to the right of the tagged particle) will asymptotically 
behave exactly like the tagged particle in Spitzer's equal mass case, cf. subsection 
above. So, the more recent and more accurate numerical results published in Boldrighini, 
Frigio, Tognetti (2002) suggesting that 

lim a 2 M =: ol « 0.74 ... (8) 

which is strictly inbetween q 2 ~ 0.627. . . and a 2 « 0.798 . . . , came as a surprise. 

The results of the present note provide substantial theoretical and independent nu- 
merical support of this surprising fact. 



4 The M —> limit of dynamics of type I 

Theorem 1. Let z G R + , u G [-1, 1], W G R be fixed, M n -> 0, V n = Mn 1/2 W n so that 
W n — * W G R, and define c := \Wz\. Choose ^ G so that for all n the dynamic 
trajectories sj:' Mn u~, z, u, V n ) and S^ i,c (uj + , u~, z) are well defined for all t G [0, oo) . 
(Mind that for any choice of z, u, W and sequences M n , V n these u^-s are of full ^ 
measure in Q^.) Then for all t G [0, oo) 

lim HSt' Mn (uj + , z, u, V„) = S^' c (u + , u~, z). 
The convergence is uniform on compact intervals of time. 

Proof. Within this proof it is convenient to describe the systems of particles as seen from 
a fixed external frame of reference: the postion of the tagged particle (in the system 
of type /) at time t is Q T t ' M , the postitions of the untagged gas particles are y±i(t), 
i = 1, 2, . . . , as given in (JSJ). 

We have to prove that in the limit described in the theorem the trajectories of the 
particles in system I converge to the coresponding trajectories in the limit system of type 
77. Note that the particles with % ^ ±1,0 follow the same dynamical rules in the two 
types of dynamics. Thus we only need to understand how the motion of the particles 
with indices ±1 can be approximated as M — > 0, if they interact with the tagged particle 
according to the rules of section 12.2.11 

As mentioned above, the particles with indices —1, and 1 have positions y~x{t) < 
Qt < 2/i (t) and velocities t>_i, V t = W t /yM and Vi(t), respectively, where v±\ and W 
are of order one. As V is very large, the tagged particle performs a full cycle: hits one of 
its neighbours, turns back, collides with the other neighbour, and gets back to its initial 
position within a very short time dt. See FigElfor insight. 

To investigate how the system evolves in this time interval dt, two successive collisions 
should be taken into account. We may assume that W > (the case of negative W is 
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t „ 



Figure 2: Successive collisions of the tagged particle of mass M <C 1 with its neighbours 

analogous), and thus the tagged particle collides first with the particle of index 1, and 
then with that of index —1. These collisions split the time interval into three smaller 
subintervals. Let us expand the formulas of (@J) in the limit as M —>■ 0. This way we may 
calculate the velocities of the three particles in the three subintervals, see Table ^ Mind 
that the order of magnitude of the incoming velocities are as follows 

V = W/VM x M" 1/2 , v±i x 1. 





1. interval 


2. interval 


3. interval 


particle 1 


Vl 


vi + 2MV + 0(M) 


vi + 2MV + 0{M) 


particle 


V 


-V + 2vi + 0(VM) 


V - 2 Vl + 2t»_! + O(VM) 


particle -1 


V-l 


V-l 


u_i - 2MV + 0(VM) 



Table 1: Velocities in the M — > approximation 



In particular, the tagged particle reverts its velocity at collisions, thus its speed is con- 
stant, more precisely, equal to \ V\ +0(1) = \ W/VM\ +0(1) throughout the investigated 
time interval. This implies dt x \[M. Furthermore, the velocities v\ and V-i remain 
0(1), and thus, the particles of index 1 and —1 remain 0(V / M)-close to their original 
positions yi and y_x in the investigated interval. In accordance with the notations of 
subsection 12+t let us introduce z — yi — y_i for brevity. By the above observations 
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the distance of the two non-tagged particles remains 0(y/M)-c\ose to z. Now we may 
calculate the leading term in dt : 



Let us denote the amount with which the velocities change at the collisions (and thus, 
during the studied time interval) by dt>i and dv-\. 
Referring to Table [T] we get 



dv 1 = 2WVM + 0(M), dv_ 1 = -2WVM + <D(M). (10) 

Formulas Q and (fTUj) altogether imply that, as M — > 0, Vi(t) and f-i(t) approach 
(piecewise) differentiable functions, and 

Vi = , V-i = . (11) 

z z 

Referring again to Table we may calculate the amount of change in the velocity of 
the tagged particle during the time interval dt. We get 



dK = 2w_ 1 -2t; 1 + C(VM), thus dW = VM(2u_i - 2v ± ) + 0(M). (12) 
Formulas and (fT^j) altogether imply 

w _ (vi-V-i)W _ zW 

z z ' 

where we have used that V\ — V-\ = yi — y~i — z. 

Integrating this differential equation we find that \Wz\ is an integral of motion. This 
means that in the M —>■ limit, c = | | is constant during a time interval when 
the tagged particle has the same neighbours, and within such an inverval, (fTTj) gives 
that the positions and velocities for the particles with indices ±1 satisfy the coupled 
differential equations 

Xl = Vi, Vi = c 2 /z 3 , 

x-i = V-i, V-i = —c 2 /z 3 . 

This is in agreement with the formula © for the potential describing the systems of type 
/J. 

Notice now that the value of c = |W / (t)z(t)| also remains constant when one of the 
neighbours of the tagged particle 'meets' another gas particle, and the neighbour is 
replaced by that new particle. At such a time moment the values of both W and z are 
unchanged. Thus for any t > |W(£)z(£)| = |W(0),z(0)| = c according to the choice of c 
in the formulation of Theorem ^ This completes the proof of Theorem ^ □ 
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Remark. Note that taking the c — > limit of the dynamics of type II, we recover the 
dynamics of type / with equal masses: the interaction between the two central particles 
becomes hard core specular collision. So, in this double limit (M — ► and than c —>■ 0)) 
the system behaves indeed as Spitzer's model, see subsection 13 . 1 1 

Recall that according to Q the measure fi 11 which is the projection of the measures 
^i,m on S £ a t e space Q 11 , is the g(c) dc-mixture of the Gibbs measures /i 7/ ' c . This 
implies that Theorem ^ has the following immediate corollary: 

Corollary. Let M — > 0. For any fixed < T < oo ; the sequence of processes [0,T] 3 
1 1 — > Q t ' converges weakly (in distribution) to the process to [0,T] 3 1 1— > Ql 1 . 

5 Numerical results on systems of type 77 
5.1 Generalities 

In this section we describe numerical investigations aimed at calculating the limiting 
variance 

a 2 := lim t -1 Var(Q t ) 

for the systems of type 71. We shall also comment on how these results are related to the 
M — ► limit of the variance for the systems of type /, as established numerically in 0j. 

These simulations of the systems of type II were done by following a number of 
particles for some fixed time T. The particles followed were those who were less than 
10T far away from the point of observation in the beginning. It is easy to check that with 
this method the probability of not following a particle that would indeed participate in 
the interaction, is negligible in all the cases we looked at. 

Numerial simulation of the dymanics of type II is relatively fast, since the equation 
of motion for the two particles interacting via the potential © can be solved explic- 
itly. (This observation is at the heart of the complete solvability of the Calogero-Moser- 
Sutherland model.) 

The simulation for time T was repeated over a sample of iV initial conditions chosen 
independently according to the appropriate stationary Gibbs distribution. From this 
sample, the empirical variance was calculated for Var(Qt) as a function of t. 

The result of a typical simulation can be seen in Figure El The solid line is the best 
linear fit for the tail, while the dashed lines have slope a 2 and a 2 , and are drawn for 
comparison. 

As we can see, Var(Q t ) does appear to be asymptotically linear. To read the limit 
lim t _1 Var(Qt) from the graph, we needn't perform such a long simulation, during 

t— >oo 

which this limit is well approached: the slope of the asymptote can be found with a 
good accuracy much sooner. Thus all the limits given in the paper are obtained using 
this technique, and the time interval for the simulation is typically between T = 10 and 
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Figure 3: Var(Qt) as a function of t in a typical simulation 

T = 50. In exchange, the size of the sample can be very big - actually, samples up to 
N = 10 7 were used. 

Finally, the statistical error of the calculated values was estimated by simply repeating 
the whole procedure about 20 times and calculating the standard deviation of the values 
obtained. 

A detailed description of the numerical simulation and the source code for the applied 
program can be found in |22j . 

5.2 The systems with fixed c 

We simulated numerically the dynamics of type II for various fixed values of the pa- 
rameter c ranging between 0.01 and 100. We started the system from samples of the 
stationary Gibbs distribution fi II,c and computed the limiting variance 

a 2 c : = lim r 1 Var(gf' c ). 

We found the c i— > a 2 c dependence of the limiting variance as shown in Figure 0] We see 
that lim c ^ °"c ~~ > which is no surprise, since in the c — > limit the system indeed 
behaves like the system of type / with M — 1, which is known to have o\ l=x = a 2 . See 
subsection 13.11 and the Remark after the proof of Theorem 1 . 

On the other hand, it is interesting to see that as c — ► oo, the limiting variance 
decreases, and even seems to approach a value near the lower limit a 2 , but not quite 
reaching this lower bound. We plan to return to this phenomenon in the forthcoming 
paper pQ. 
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Figure 4: c- dependence of the limiting variance for systems of type II 



5.3 The mixed system 

We computed the numerical value of the limiting variance for the mixture of dynamics 
in two different ways. 

First, we computed numerically the value of 

a 2 ^ x := lim r^Var (Qf). 

We have done it in the following way: we sampled the initial conditions u 11 according 
to the distribution /i 7/ and, independently, a standard normal variable W. Then we 
computed c := \Wz\, where z was the distance between the two central particles in 
the initial configuration uj n . This random value c served as strength parameter in the 
interaction potential ©, with which the dynamics <S/ /,C was computed. 

Second, using the data obtained for of in the fixed c computations (see subsection l5.2|) . 
we computed the mixture 

poo 

^mix,2 : = / °cP(c) dc, 
JO 

which, of course, in principle must give the same value as the previous computation. 
Indeed, in the two cases we obtained the numerical values 

a^ ix l = 0.736 ± 0.003, cx^ ix 2 = 0.740 ± 0.003. 

This result is very interesting, since it coincides exactly (well within statistical error) 
with the M — > limit of the variance calculated numerically in see (jH}- This means 
that there is indeed continuity in the limiting variance as M — > 0. 
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We remark that the function c h- > o 2 c shown in Figure 0] can be fit with amazing 
accuracy by the function of the simple form a 2 = ^rry + ^2, where A\ = 0.796, 

A 2 = 0.638, c = 1.981, p = 0.792. 
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